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Abstract 

More than a dozen binary systems are now established as sources of variable, high energy (HE, 0.1 — 100 GeV) 
gamma rays. Five are also established sources of very high energy (VHE, ^ 100 GeV) gamma rays. The mechanisms 
behind gamma-ray emission in binaries are very diverse. My current understanding is that they divide up into 
four types of systems: gamma-ray binaries, powered by pulsar rotation; microquasars, powered by accretion onto 
a black hole or neutron star; novae, powered by thermonuclear runaway on a white dwarf; colliding wind binaries, 
powered by stellar winds from massive stars. Some of these types had long been suspected to emit gamma rays 
(microquasars), others have taken the community by surprise (novae). My purpose here is to provide a brief 
review of the current status of gamma-ray emission from binaries, in the context of related objects where similar 
mechanisms are at work (pulsar wind nebulae, active galactic nuclei, supernova remnants). 

Resume 

Emission gamma des systemes binaires. Plus d’une douzaine de systemes binaires sont maintenant identihes 
comme des sources variables de rayonnement gamma de haute energie (0.1 — 100 GeV). Cinq systemes binaires 
sont egalement des sources de rayonnement gamma de tres haute energie (J? 100GeV). Les processus menant a 
remission de ce rayonnement sont varies. On pent neanmoins diviser ces systemes en quatre grandes classes : les 
binaires gamma, dont le moteur est la rotation d’un pulsar; les microquasars, dont le moteur est I’accretion sur 
un trou noir on une etoile a neutrons; les novae, dont I’energie provient de la combustion nucleaire a la surface 
d’une naine blanche; les binaires a collision de vents, qui dissipent I’energie cinetique de vents d’etoiles massives. 
On soupQonnait depuis longtemps que certaines classes de systemes binaires devaient etre associees a de remission 
gamma (les microquasars); pour d’autres classes, la decouverte d’emission gamma a ete une surprise (les novae). 
Je propose ici une breve revue de nos connaissances sur remission gamma de haute et tres haute energie dans 
les systemes binaires, avec le souci de souligner les liens avec d’autres objets astrophysiques dans lesquels des 
processus similaires sont a I’oeuvre (nebuleuses de pulsar, noyaux actifs de galaxie, restes de supernova). 
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1. Introduction 


Detections of binaries punctuated high energy and very high energy gamma-ray astronomy in the 1970s and 
1980s. Binaries had been discovered to be the brightest sources in the newly-opened X-ray sky and were thus 
likely candidates for emission at even higher energies. Most of these detections did not withstand the test of time 
as independent confirmation could not be obtained [1]. Yet, they fostered great interest and participated in the 
development of the field from the 1980s onwards. 

The latest generation of gamma-ray instruments, from the mid-2000 onwards [2, 3], have brought back binaries 
to the forefront. Many binaries are now established sources of high energy (HE, 0.1 — 100 GeV) and very high 
energy (VHE, 100 GeV) gamma rays: I have listed in Table 1 the names and main characteristics of those that 
are reported as secure gamma-ray detections and display variability on a timescale that can be related to the 
binary nature of the source. I have also listed one or two references to the most recent HE and/or VHE studies as 
a point of entry into the litterature. Some of these binary gamma-ray sources are of the type that had long been 
suspected to emit gamma rays (microquasars), others have taken the community by surprise (novae). There is no 
simple extrapolation between the HE and VHE domains. At present, the only binaries detected in VHE gamma 
rays appear to be composed of a pulsar and a massive star. There is much more diversity in HE gamma rays, 
perhaps because the higher fluxes and the instrumental sensitivity at these energies give access to more objects. 

My purpose here is to sketch a status of gamma-ray emission from binaries as of early 2015 and how it is 
understood in the broader context of related objects and high-energy emission models. More in-depth accounts 
of the observations and models may be found in other reviews [4, 5, 6, 7]. My current understanding is that the 
binaries detected in HE or VHE gamma rays divide up into four classes, illustrated in Eigure 1: 

- gamma-ray binaries, powered by pulsar rotation; 

- microquasars, powered by accretion onto a black hole or neutron star; 

- novae, powered by thermonuclear runaway on a white dwarf; 

- colliding wind binaries, powered by stellar winds from massive stars. 

The last three classes were well-established prior to the detection of gamma rays. Gamma-ray binaries are still 
rather new and, although some of the systems had been known for some time, the motivation to treat them as 
a separate class is an outcome of the last decade of gamma-ray observations. The contours of this class are not 
fixed: sometimes the term “gamma-ray binary” is used in the literature to designate all binaries that have been 
detected in HE or VHE gamma rays (all the above), sometimes it encompasses only the systems listed as high- 
mass gamma-ray binaries in Table 1, sometimes it adds the microquasars to the high-mass gamma-ray binaries. 
I think the first definition is too broad, the second turns out to be too restrictive, the third mixes objects with 
very different phenomenology and underlying powering mechanisms. 

A distinguishing feature of a gamma-ray binary is that emission above 1 MeV dominates its spectral energy 
distribution (see example in Figure 2), disregarding the black body like component from the companion star [6]. 
Another distinguishing feature of the gamma-ray binaries listed in Table 1 is that they all have variable gamma- 
ray emission, sometimes modulated on the orbital period. This definition encompasses both binaries in which the 
companion to the compact object is a high-mass (^ 10 M©) star and those where the companion is low-mass 
(^ IM©). The ensuing subdivision into low-mass and high-mass gamma-ray binaries, like X-ray binaries, actually 
makes sense as I shall explain in section 2. Radio pulsations identify the compact object as a neutron star in many 
of these systems. Theoretically, the prevailing idea is that gamma-ray emission from these systems is powered by 
the rotation of the pulsar magnetosphere. 

X-ray binaries are powered by accretion of matter from a normal star onto a stellar-mass black hole or a 
neutron star. Infalling matter forms an accretion disk around the compact object, where its gravitational energy 
is gradually released. In some X-ray binaries, part of this energy ends up powering a jet i.e. a collimated ejecta 
with relativistic speeds, similar to those seen in active galactic nuclei (radio quasars, blazars), hence the term 
“microquasar” [8]. Most of the radiated energy is released in X-rays. The prevailing idea is that the gamma-ray 
emission is associated with the formation of the relativistic jet. 

Novae events are due to runaway thermonuclear burning of material deposited on the surface of a white dwarf. 
Novae occur in cataclysmic variables, binaries with a white dwarf accreting from a stellar companion [9]. Depending 
upon the mass transfer rate, it can take tens to thousands of years for the accumulated material to reach the critical 
pressure that triggers thermonuclear burning. The prevailing idea is that the gamma-ray emission is associated 
with the ejection of part of the material when runaway burning starts. The rest of the material inflates into a 
large enveloppe around the white dwarf until, after days to weeks, nuclear burning ceases and the pressure is 
insufficient to maintain the enveloppe. Most of the radiated energy is released in the optical band where they are 
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name 

binary components 

Porb (d) 

HE 

VHE 

refs (*) 

notes 


(high-mass) gamma-ray binaries 

PSR B1259-63 

pulsar 

Be 

1236.7 



[12, 13] 

47.7 ms 

HESS J0632-I-057 

? 

Be 

315 



[14. 15] 


LS I +61°303 

? 

Be 

26.5 

■/ 


[16. 17] 

magnetar ? 

IFGL J1018.6-5856 

? 

o 

16.6 



[18, 19] 


LS 5039 

? 

o 

3.9 

■/ 


[20. 21] 


(low-mass) gamma-ray binaries (f) 

XSS J12270-4859 

pulsar 

red dwarf 

0.29 



]22, 23] 

1.7 ms 

PSR J1023+0038 

pulsar 

red dwarf 

0.20 

■/ 


[24] 

1.7 ms 

2FGL J0523.3-2530 

? 

red dwarf 

0.69 

■/ 


]25. 26] 


PSR B1957+20 

pulsar 

brown dwarf 0.38 



[27] 

1.6 ms 

PSR J0610-2100 

pulsar 

brown dwarf 0.29 



[28] 

3.8 ms 

PSR J1311-3430 

pulsar 

brown dwarf 0.065 

■/ 


[29, 30] 

2.6 ms 


microquasars (X-ray binaries) 

Cyg X-3 

black hole ? 

Wolf-Rayet 

0.20 



[31, 32] 


Gyg X-1 

black hole 

o 

5.60 

■/ 

? 

]33, 34] 



novae 


V407 Cyg 

white dwarf 

red giant 

14000 ? 


]35, 36] 

N Cyg 2010 

V1324 Sco 

white dwarf 

red dwarf 

0.07 ? 


[37] 

N Sco 2012 

V959 Mon 

white dwarf 

red dwarf 

0.30 


[37] 

N Mon 2012 

V339 Del 

white dwarf 

red dwarf 

0.13 ? 


[37] 

N Del 2013 

V1369 Gen 

white dwarf 

red dwarf 

? 

■/ 

[38] 

N Cen 2013 

colliding wind binary 

Eta Car 

LBV 

O/WR ? 

2014 


[39, 40] 



* I only give one or two recent references as entry points to the HE/VHE litterature. 
t not including another > 50 Eermi-LAT pulsars in binaries. 

Table 1 

Binary systems identified with sources of variable HE or VHE gamma-ray emission. 

sometimes visible to the naked eye. 

Colliding wind binaries are composed of massive O or Wolf-Rayet stars. These stars have masses in excess of 
20 Mq and large luminosities 10'* —10® Lq . The radiation pressure drives strong winds from the stars with terminal 
velocities « 1000 — 2000 km s“® and mass loss rates of 10“® — lO“®M 0 yr“^. The collision between the stellar 
winds in a binary forms a shock structure where the released kinetic energy heats the gas to ~ 10^ K, which emits 
principally X-rays [10, 11]. Some of the energy may be channeled into accelerating particles at the shock, as in 
supernovae remnants, leading to gamma-ray emission. 

The detection of gamma rays has brought new light, literally, to the mechanisms in operation in these systems. 
The mechanisms behind gamma-ray emission in binaries are very diverse (see Figure 1 for illustrations of these). 
In novae, as will be explained in section 4.2, the physics has close analogies to that of supernova remnants; the 
gamma-ray emission from microquasars involves relativistic jets as active galactic nuclei and gamma-ray bursts 
do; many binaries host pulsars whose rotating magnetic field ends up powering strong non-thermal emission. 
Table 1 makes it clear that the binaries emitting gamma rays do not constitute a class of objects united by a 
common astrophysical model, as for instance supernova remnants, but by a common phenomenon: orbital motion. 
Orbital motion constraints the size and type of the binary components, and allows measurements to probe different 
conditions or vantage points. The legitimacy of theories developed for gamma-ray sources in other astrophysical 
contexts can thus be put to the test in the peculiar environment of binaries. 

In the following, I have tried to emphasize how gamma-ray observations of binaries provide such tests. Thus, 
section 2 sets gamma-ray binaries in the context of pulsars, section 3 relates microquasars to the physics of 
relativistic jets, and section 4 examines the gamma-ray detections of novae and colliding winds binaries in the 
context of particle accretion at non-relativistlc shocks. 
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Figure 1. Sketch of the types of gamma-ray emitting binaries, (a) A high-mass gamma-ray binary, with a rotation-powered puisar. 
The filled in region represents the shocked pulsar wind and the shocked stellar wind (panel adapted from [41]). (b) A microquasar, 
with its accretion disk -h corona, and the relativistic jet. (c) A nova, with the expanding ejecta and enveloppe around the white 
dwarf, (d) A colliding wind binary, with two stellar winds interacting. The sketches are not to scale. 


2. Gamma-ray binaries and pulsars 
2.1. PSR B1259-63 as a binary PWN 

The detection of PSR B1259-63 with the H.E.S.S. arrays of Cherenkov telescopes was the first detection of a 
binary with the latest generation of gamma-ray instrumentation [42]. As the name implies, PSR B1259-63 is a 
binary system harbouring a pulsar. Its 48.8 ms pulsations were detected with the Parkes radio telescope whilst 
surveying the sky for new pulsars [43]. The radio pulsations can be tracked very accurately, allowing for a precise 
determination of the orbit and the measurement of the rate at which the pulsar rotation slows down. The slow down 
implies rotational kinetic energy is lost at a rate E ~ 8x 10^® ergs“^ (8 x 10^® W), which is a high spindown power 
for a pulsar. This energy loss basically corresponds to the power required to maintain the magnetic field in rotation 
with the neutron star. At zeroth order the pulsar behaves like a rotating dipole, enabling the pulsar spindown 
power to be related to the magnetic field. In PSR B1259-63 the magnetic field of the pulsar is B ~ 3.3 x 10^^ G 
(3.3 X 10^ T). These characteristics make PSR B1259-63 a rather typical rotation-powered young pulsar [44[. 

The VHE and HE gamma-ray emission are not pulsed but are only detected close to periastron passage, when 
the pulsar is closest to its companion. The companion is a massive star with a mass « 30 Mq [45]. In addition 
to the fast, radiatively-driven stellar wind, spectroscopy shows emission lines arising from a slowly outflowing 
disk of material in Keplerian rotation around the equator of the star. The origin of these circumstellar disks, 
characteristic of “Be stars”, is not known but thought to be due to the star rotating close to breakup velocity 
[46]. The 3.4 year orbit has an eccentricity of 0.87 that takes the pulsar from 13 AU to 0.9 AU from the star. At 
periastron, the pulsar is about 20 stellar radii away from the star. The density of material in the circumstellar 
disk and stellar wind is then high enough to eclipse the radio pulsations for « 40 days around periastron due to 
free-free absorption [47]. 

In rotation-powered pulsars, most of the energy lost by pulsar spindown powers a tenuous, highly relativistic 
wind carrying e^e“ pairs and electromagnetic energy away beyond the pulsar light cylinder radius Rl [48, 49], the 
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radius at which the rotation velocity of material tied to the pulsar magnetic field lines reaches the speed of light 
and must then decouple {Rl = cP/2t: where P is the pulsar spin period). When the pulsar is isolated, the pulsar 
wind propagates freely until it is contained by the interstellar medium or the ejecta from its parent supernova, 
forming a shock structure. These are called pulsar wind nebulae (PWN) [50, 51]. In PSR B1259-63, the pulsar 
wind propagates in the much denser environment provided by the circumstellar disk and stellar wind. As a result, 
the pulsar wind is (at least partly) confined, forming a shock structure with the surrounding stellar material on 
a much smaller scale than in a PWN [52]. 

The distance to which the pulsar wind propagates freely is estimated by equating the pulsar wind ram pressure 
to the pressure of the surrounding material. If the surrounding (ram) pressure is due to an isotropic, coasting 
stellar wind, the distance R to the pulsar is 


R 1 MwUw 

— ~ -TTTT With ri = — ^- 

d 1 + 171/2 


( 1 ) 


where d is the orbital separation, E is the pulsar wind power, and the stellar wind mass loss rate and 
velocity (see [ 6 ] for details and references). A stellar wind with a strong ‘thrust’ M„v„ contains the pulsar wind 
close to the neutron star (77 ^ 1 +> R/d <C 1). Inversely, a strong pulsar wind pushes back the stellar wind all the 
way to the surface of the star (17 <C 1 => R/d ~ 1). 

Here, as in PWN [51], the shock is thought to channel a fraction of the pulsar wind energy to particles that then 
radiate gamma rays. The process may involve randomization of the particles initially frozen with the magnetic field 
in the pulsar wind, diffuse shock acceleration, or shock-driven reconnection. The latter is increasingly scrutinized 
since it has become clear that highly relativistic pair plasmas are poor at shock-accelerating electrons to very high 
energies unless they are very weakly magnetized (probably too much so to be reconciled with observations and 
with the requirement that the pulsar wind starts at the light cylinder by carrying a lot of magnetic energy) [53]. 
Anyhow, we know that the process must be efficient because a large fraction of the spindown power ends up as 
non-thermal radiation: the peak luminosity in HE gamma rays from PSR B1259-63 is nearly equal to its spindown 
power [ 12 ] ! 

Besides PSR B1259-63, there are four other known pulsars with massive star companions: PSR J1740-3052, PSR 
J1638-4725, PSR J0045-7319, and PSR J2032+4127. The first three are probably too faint in gamma rays to be 
detectable because of their lower spindown power and larger distance than PSR B1259-63 [ 6 ]. The last one, PSR 
J2032+4127, has 143 ms pulsations detected in radio and HE gamma rays and has been tentatively associated with 
VHE source TeV J2032+4130 [54]. PSR J2032+4127 is actually in a long (20-30 year), highly eccentric (e « 0.95) 
orbit with a 15 Mq star [55[. Like PSR B1259-63, variable gamma-ray emission may be detectable in early 2018 
when periastron passage will bring the pulsar within 2 AU (3 x 10^^ m) of its companion star. 

Four other binaries, listed in Table 1, currently have observational properties that are similar to PSR B1259-63 
(modulated HE or VHE gamma-ray emission, massive star companion, etc). Based on this, the binary PWN 
scenario described above for PSR B1259-63 has been widely adopted for those four systems, even if some contest 
this interpretation, favoring a microquasar scenario [56]. The basic issue is that it has not been possible to detect 
radio pulsations, most likely because these binaries have much shorter orbital periods than PSR B1259-63. The 
pulsar is so deeply embedded within the dense stellar wind that the radio signal cannot propagate out [57]. Still, 
two magnetar bursts in X-rays have been detected from the direction of LS I +6T‘’303, within « 2 arc minutes of 
the source, strengthening the case for a pulsar in this system [58]. I am of the opinion that the array of evidence 
strongly favours the binary PWN scenario and refer to [ 6 ] for a detailed discussion. 


2.2. Low-mass gamma-ray binaries 

Gamma-ray emitting pulsars are also found in association with low-mass companions: red dwarfs, brown dwarfs, 
white dwarfs, other pulsars. For instance, HE gamma rays were detected from the famous double pulsar PSR J0737- 
3039A/B [59]. As of November 2014, 55 of the 161 pulsars in the Fermi-LAT catalog of gamma-ray pulsators ^ 
are actually in binaries, far outnumbering the systems listed in Table 1. These pulsars have millisecond pulse 
periods and belong to the class of old pulsars that have been “recycled" i.e. spun back up by accretion of material 
from the companion onto the neutron star [44]. In nearly all of these systems, the only difference with an isolated 
gamma-ray pulsar is that the pulse frequency is modulated by the orbital motion of the pulsar [60]. 


1. https://confluence.slac.stanford.edu/display/GLAMCOG/Public+List+of+LAT-Detected+Gamma-Ray+Pulsars 
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However, there are a few of those systems where their presence in a binary has an impact on the gamma- 
ray emission. These are the five systems that I have listed in Table 1 as low-mass gamma-ray binaries, because 
gamma-ray emission dominates their spectral energy distribution, they show variability linked to their binary 
nature (albeit often tentatively), and they also involve interacting pulsars. 

The HE gamma-ray emission from PSR B1957+20, PSR J0610-2100, PSR J1311-3430, and 2FGL J0523.3-2530 
shows some evidence for a modulation on the orbital period [27, 28, 26]. PSR B1957-I-20 is the original “black 
widow” pulsar, named as such because the pulsar wind in this system impinges directly onto its very low mass 
companion 0.1 M©), ablating material from the star (see [61] for a review). There is no strong enough stellar 
wind or accretion flow present to stop the pulsar wind {rj 1). Part of the gamma-ray spectrum is not variable and 
attributed to pulsar magnetospheric emission, as in isolated pulsars, while the emission above « 2 GeV is variable 
and attributed to the pulsar wind. The processes are identical to those invoked for the “high-mass” gamma-ray 
binaries {e.g. [62, 63[). The modulation is detected in only part of the Fermi-LAT data on 2FGL J0523.3-2530, 
raising the interesting possibility that orbital modulations may also be sporadically present in some of the other 
50 or so pulsars in binaries mentioned above. 

PSR J1023+0038 and XSS J12270-4859 are nearly twin systems (companion, orbital period, pulse period), part 
of the “redback” class of pulsars (basically black widows with companion masses of a few tenths instead of a few 
hundredths of solar mass, see [61]). Both are intermittent HE gamma-ray sources and both show intermittent 
signs of accretion, manifest as an increase in X-ray and optical/UV flux, stronger variability, and double-peaked 
optical emission lines typical of accretion disks. Both are sometimes detected as radio pulsars, and sometimes 
detected as X-ray pulsars [64, 65]. The neutron star seems to transition back and forth between accretion-powered 
and rotation-powered regime, even combining both at times ! Another such “transitional millisecond pulsar” is 
PSR J1824-2452I in the globular cluster M28 [66]. Although gamma-ray emission is detected from M28, most of 
the flux seems to arise from another pulsar in the globular cluster [67] so I have not included it in Table 1. The 
Fermi-hAT sources 3FGL J1544.6-1125 and IFGL J1417.7-4407 also very likely belong to this sub-class of sources 
[68, 69]. I have not included them in Table 1 because gamma-ray variability (pulsations or flux changes) has yet 
to secure the identifications. 

These transitions are a key verification that old, millisecond radio pulsars have been spun up by accretion. 
They also represent a challenge since the models predicted that the pulsar wind, once the radio pulsar turned on, 
would quench accretion and forbid a return to an accreting state [70]. Curiously, in PSR J1023+0038, the HE 
gamma-ray flux increased by a factor 5 when the source entered an active X-ray state, whereas it was expected 
that the gamma-ray emission associated with the radio pulsar would disappear with the onset of accretion. A 
possibility is that the pulsar mechanism is not turned off and that the additional light from the active accretion 
disk, which does not extend to the neutron star, leads to higher levels inverse Compton gamma-ray emission from 
high-energy particles in the pulsar wind [71]. 


2.3. Gamma-ray binaries and the physics of pulsar winds 

Binaries give access to pulsar winds on smaller scales than in PWN since the termination shock distance is a 
fraction of an AU (the orbital separation), compared to a fraction of a parsec in a PWN like the Crab nebula. In 
units of the light cylinder radius Rlc (the natural scale for the pulsar wind), both low and high-mass gamma-ray 
binaries probe scales of ~ IO'^Rlc instead of ~ 10®i?LC for an isolated pulsar [72[. The star provides a source 
of photons that can scatter off the high energy electrons in the pulsar wind and the shocked region. The stellar 
radiation field is very well known so the resulting up-scattered gamma-ray spectrum constrains the location and 
properties of the high energy electrons. However, this is complicated by the anisotropy of the radiation processes 
involved (inverse Compton emission, pair production) and by relativistic Doppler boosting (the shocked wind 
maintains mildly relativistic speeds) which, together with the varying orbital separation, contribute to modulate 
the non-thermal emission up to the gamma-ray regime [6[. 

Uncertainties in the geometry of the shock region also complicate the interpretation. Estimating the exact 
location of the termination shock is sometimes difficult because the interaction produces a complex, 3D, time- 
dependent morphology, notably when a circumstellar disk is present. This has begun to be addressed through 
numerical simulations [73, 74, 75, 76, 77]. Studying systems where the massive star is an O star, without a 
circumstellar disk, simplifies the geometry and may be more amenable to interpretation. Yet, the observations 
pose an increasing number of challenges to models. 

First, the gamma-ray spectra show unexpected complexity with two, if not three, distinct emission components 
in high-mass gamma-ray binaries. For example, in LS 5039 (see Figure 2), the MeV spectrum (CGRO-COMPTEL) 
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does not connect to the GeV spectrum (Fermi-LAT), itself clearly distinct from the TeV spectrum (H.E.S.S.) [78]. 
The gamma-ray fluxes are modulated on the orbital period, but with different phasing (GeV vs TeV and MeV, see 
Figure 2). Cross-calibration issues between instruments are not suspected, so this signals different populations of 
high-energy particles and/or emission mechanisms.Figure 2 shows a simple one-zone model where electrons with 
a E~^ distribution in energy radiate synchrotron and upscatter blackbody photons from the star (at 0.2 AU), 
illustrating why the GeV emission requires an additional population of electrons ^ . Various possibilities have been 
proposed for the origin of the GeV emission. An exciting proposition is that it is the inverse Compton emission 
of pairs in the pulsar wind [80, 81, 82, 41[. The pulsar wind is usually ‘dark’ in PWN because the particles are 
frozen in with the magnetic field and the background radiation field is weak. In binaries, the companion behaves 
as a torch illuminating the pulsar wind, instilling hope of a direct access to its energetic and particle content. 

Second, the gamma-ray lightcurves show much more complex behaviour than straightforward ‘geometrical’ 
radiative models usually account for. These models are based on the change with orbital phase of the binary 
geometry seen by the observer. Typically, strong gamma-ray emission will be expected when the high-energy 
electrons are close to the source of seed photons and scatter back the radiation in the direction of the observer 
{e.g. [83, 84[). For instance, in PSR B1259-63, the gamma-ray emission due to inverse Gompton scattering should 
peak shortly before periastron passage if most of the high-energy electrons reside close to the pulsar and most 
of the seed photons come from the companion star. However, the HE lightcurve actually peaks a month after 
periastron passage, when the scattering geometry is unfavourable [41[. The HE luminosity also implies a very 
high efficiency in converting the spindown power into gamma rays that turns out to be hard to achieve, unless 
some relativistic boosting is involved [85[. Another difficulty raised by the gamma-ray lightcurves are the changes 
in the HE and VHE modulation of LS I +61‘’303 on timescales both longer and shorter than the orbital period 
[17, 86[. Modelling will be difficult if a lot of the gamma-ray variability happens to be related to fluctuations 
in the circumstellar environment of the pulsar (clumps in the stellar wind, cycles in stellar mass-loss rate, spiral 
waves in the Be disk, etc [7, 87]). Observations on a long time base are required to extract the average behaviour, 
and thus identify the regular variations due e.g. to the periodic compression of the pulsar wind (when the orbit 
is eccentric) from the more irregular ones due to the companion’s weather. 

Finally, robust conclusions on pulsar winds require good knowledge of the binary system parameters: orbital 
period, eccentricity, orientation, but also stellar wind mass loss rate, velocity and, ideally, the pulsar spin period 
and spindown power [88[. This represents a major observational effort, and not all of the parameters are easily 
accessible in all systems. In light of this, increasing the number of known gamma-ray binaries will help identify 
the more accessible systems, in addition to contributing to our knowledge of this population. 


3. Microquasars 

Microquasars, as scaled-down versions of active galactic nuclei (AGNs), have long been on the list of likely HE 
and VHE gamma-ray emitters. Relativistic jets launched from the close environment of black holes in AGNs are 
known to produce powerful, variable HE and VHE gamma-ray emission in blazars [93[. Furthermore, the detection 
in several microquasar jets of synchrotron emission extending from radio to X-ray energies unambiguously proves 
that TeV electrons are present on parsec scales, as in AGN jets [94]. 

There are currently two microquasars detected in HE gamma rays with both AGILE and the Fermi-LAT. 
Cyg X-3 [95, 96, 31[ and Cyg X-1 [97, 98, 33[ are well-known systems composed of a black hole candidate in orbit 
with a massive star. Their X-ray luminosity is around 10®® to 10®^ W, 3 to 4 orders-of-magnitude greater than in 
gamma-ray binaries with massive companions. This emission arises from the accretion disk and its corona, with a 
contribution from the jet in some models [99[. The presence of relativistic jets in both systems is well-documented 
from radio observations, tracing non-thermal synchrotron emission from the ejecta. The activity of these jets is 
closely related to the X-ray spectral state [8[. The level of gamma-ray emission does not exceed « 1% to 10% 
of the X-ray luminosity. VHE observations have yielded only upper limits, except for a single flare event during 
MAGIC observations of Cyg X-1 with a significance of 4.9a (pre-trial) [34[. Subsequent long campaigns by the 


2. The magnetic field intensity {B ~ 0.04 T) is set by the ratio of synchrotron to Compton fluxes [53]. The gamma-ray spectrum 
fails to reproduce the hard H.E.S.S. spectrum, mostly because it is calculated in the isotropic approximation whereas it should take 
into account that the incoming photons from the star come from a specific direction as well as relativistic effets [79]. The slope 
of the synchrotron emission fits well with the X-ray to MeV spectrum, but note that this is not self-consistent since very strong 
cooling is expected at the high-energy end. 
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Figure 2. Spectral energy distribution (left) and orbital modulation (right) of the 0.1-10 GeV (blue squares, [89]) and TeV (red 
circles, [20]) flux from the gamma-ray binary LS 5039. The spectrum includes data from Suzaku [90], BATSE INTEGRAL [91, 92], 
GOMPTEL [78], Fermi-LAT [21], H.E.S.S. [20], in order of increasing frequency. A simple one-zone model is adjusted to the data 
for illustration purposes, including synchrotron and inverse Compton emission from a single population of electrons (see section 2.3 
for details). 


MAGIC and VERITAS collaborations did not detect another such VHE flare (hence the question mark in Table 
1) [100, 101]. Cyg X-3 has not been detected in VHE gamma rays [102]. 

The major feature of the HE gamma-ray emission is that it is related to the spectral state of the binary. 
Gamma rays are detected only when radio emission is present, clearly associating this high-energy emission to 
the relativistic jet. Spectral modeling also supports this view: the HE gamma-ray emission is distinct from the 
non-thermal power law spectral tails detected up to several MeV in both systems. The GeV emission is not a 
simple extrapolation of the MeV emission. The MeV tails probably arise from electrons in the corona while the 
HE gamma-ray emission arises further out ]103, 104]. The orbital modulation of the HE emission from Cyg X-3 
also constrains the location of the particles in the jet to be at a distance comparable to the orbital separation, 
very far from where most of the gravitational energy is released [105, 106] ! This interpretation assumes that the 
modulation arises from inverse Compton scattering of photons from the massive star by electrons in the jet. HE 
gamma ray emission too close to the black hole would be strongly attenuated by pair production with the X-ray 
photons from the accretion disk [107, 108]; HE gamma ray emission too far from the black hole would hardly be 
modulated and would require too high a power in electrons to generate the gamma-ray flux against the decreased 
stellar radiation density ]105]. 

Particle acceleration at a recollimation shock within the jet has been proposed as an explanation for the distant 
gamma-ray emission site in Cyg X-3. Such shocks occur when the pressure from the expanding jet has dropped to 
match the pressure of the surrounding medium. Recollimation shocks may account for knots of emission in AGN 
jets [93], notably in M87 ]109], as particles in the jet are re-energized when they pass through such a shock. The 
Wolf-Rayet companion in Cyg X-3 has a mass loss rate and the tight orbit (4.8 hr) places the compact object 
only one stellar radius or so away from the star. The enormous ram pressure from the stellar wind is sufficient to 
cause the jet to recollimate within the binary system in Cyg X-3, though not necessarily in Cyg X-1 (which has an 
O type companion). Attributing the gamma-ray emission to a specific location may not be necessary: in models 
where non-thermal electrons are continuously accelerated and cooled by inverse Compton scattering on stellar 
photons, the jet gamma-ray emission naturally peaks at a height comparable to the orbital separation ]110]. 

Despite extensive searches, HE or VHE gamma-ray emission has not been detected yet from other X-ray binaries 
]111, 100, 101, 112, 113], including systems with low-mass companions ]114, 115, 116] and accreting neutron stars 
with high-mass companions ]32, 117]. In the latter, there is no radio emission from a relativistic jet and no pulsar 
wind. In the former, despite the presence of relativistic jets attested by radio observations, the radiation field from 
the low-mass star is much weaker. Particles in the jet may still be able to upscatter their synchrotron radiation 
to gamma-ray energies, but with a lower efficiency and closer to the base of the jet, where pair production may 
attenuate most of the emission. At very large distances, the dissipation of the jet kinetic energy in the interstellar 
medium though a shock may lead to detectable gamma-ray emission. An example is the recent report of HE 
gamma-ray emission from the microquasar SS 433 ]118]. This detection, in a region of complex diffuse gamma-ray 
emission, remains to be confirmed. 
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Figure 3. Multiwavelength lightcurve of the microquasar Cyg X-3 from 2010 September 21 to 2011 June 8, around the time of a 
giant radio flare (grey band shows the onset). The panels show (a) hard X-rays from 5'if;^/f/BAT, (b) soft X-rays from RXTE/ASM. 
and MAXI, (c) Fermi-hAT HE gamma rays, (d) 15 GHz radio flux (note the log scale). Roughly, the soft X-rays trace the accretion 
disk, hard X-rays trace the disk corona, gamma rays trace particle acceleration in the jet, and the radio trace the non-thermal 
emission from the relativistic jet. Figure reprinted from [31], with permission of Oxford University Press, on behalf of the RAS. 


Energy released by accretion also powers relativistic jets and gamma-ray emission in AGN and in gamma-ray 
bursts [93]. Yet, it is only in microquasars that the release of energy in the different channels can be followed 
simultaneously, on reasonable timescales (changes occur on timescales of hours to days), albeit at the price of 
complex multiwavelength campaigns combining multiple facilities [31]: soft X-ray trace the accretion disk, hard 
X-rays trace its corona, radio waves trace the large-scale relativistic jet, and gamma rays trace the non-thermal 
particle injection during ejection (Figure 3). Although rare, gamma-ray emission from microquasars thus provides 
unique insights into the connection between accretion and ejection. 


4. Novae and colliding wind binaries 

4.1. Eta Car and the missing gamma-ray emission from colliding wind binaries 

Colliding wind binaries (CWB) involve Interacting winds, exactly like (high-mass) gamma-ray binaries (Fig. 1) 
except that the winds both originate from massive stars, hence the shock is non-relativistic [10]. There is actually 
an evolutionary sequence from colliding wind binary, to gamma-ray binary (once one of the massive stars explodes 
and forms a neutron star), to X-ray binary (once the neutron star rotation has slowed down and the pulsar wind 
pressure is insufficient to prevent accretion). Non-thermal radio synchrotron emission has established that particle 
acceleration occurs in CWB [119]. The large photon and matter densities in these systems should easily lead to 
gamma-ray emission from inverse Compton scattering or pion-production [120]. The gamma-ray emission would 
then inform on the maximum acceleration energy, the magnetic field, and the fraction of the wind kinetic energy 
converted to non-thermal energy at the shocks. 

Despite expectations, gamma-ray emission is currently reported in only one CWB: Eta Car [121, 122, 40, 39]. 
Eta Car is composed of two massive stars, one a gigantic « 100 Mq luminous blue variable and the other an O 
or Wolf-Rayet type star, in a 5.5 year orbit [123[. HE gamma rays are detected using AGILE and Fermi-hAT up 
to 100 GeV, though the system has not yet been detected at VHE by ground-based Cherenkov arrays, suggesting 
a spectral cutoff around a few 100 GeV [40]. The gamma-ray spectrum consists of two components: soft emission 
cutting off below 10 GeV and a hard power-law tail above 10 GeV, modulated on the orbital period [39]. 
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The gamma-ray luminosity of Eta Car represents only « 0.2% of the kinetic power available from the stellar 
winds (« 6 X 10^° W). Curiously, no gamma rays have been detected from other nearby CWB with comparable or 
even higher kinetic powers [124], so their gamma-ray efficiency must be lower. Perhaps this is due to circumstances 
in the binary environment that reduce the high-energy emission {e.g. strong radiative losses, or the shocked fraction 
of the wind is small). Perhaps this is signaling that diffusive acceleration is impeded in CWB compared to supernova 
remnants {e.g. because the CWB shock structure is so unstable that the particles do not see the sharp interface 
required for Fermi acceleration). 


4.2. Novae 

In the initial stages of a nova eruption, a fraction of the white dwarf envelope is ejected while the rest, supported 
by the pressure from thermonuclear burning at the base, inflates to ~ 100 R© on a timescale of a day, followed 
by a phase of strong radiatively-driven mass loss lasting tens of days [125]. Novae are powered by nuclear burning 
so low-energy (~1-10 MeV) gamma-ray emission from radioactive decay products has long been expected from 
these systems. However, gamma-ray emission above 100 MeV came as a surprise, because there was no compelling 
reason - observational or theoretical - to expect a lot of particle acceleration to occur (but see ]126]). 

Five novae have now been detected in HE gamma rays with the Fermi-LAT ]37]. The hrst, V407 Cyg, occurred 
in a binary composed of a red giant in a long orbit around the white dwarf (such systems are called “symbiotic” 
binaries) ]35]. Red giants lose signihcant amounts of matter, ~ 10~^ to 10“® M© yr“^, in slow stellar winds with 
velocities of ~ 20kms“^. Thus, the nova ejecta had to propagate in a dense circumstellar environment, forming 
an expanding shock on scales comparable to the orbital separation (10 AU), where particles can be accelerated 
to high energies. The system behaves as a small-scale analogue of a supernova remnant ]127]. Another symbiotic 
nova (not listed in Table 1), V745 Sco, yielded only a weak 2 to 3cr signal in February 2014, probably because 
of its 7.8 kpc distance, to compare to 2.7 kpc in V407 Cyg. Many symbiotics are recurrent novae and we can 
look forward to e.g. the eruptions of RS Oph (4.2 kpc) in 2021 (H.E.S.S. observations were attempted a couple of 
months after the last eruption in 2006, probably too late) and of T Cr b in 2026 (0.9 kpc) [128] ! 

After the detection of V407 Cyg, the presence of a red giant wind seemed to be the necessary condition for shock 
formation and gamma-ray emission in a nova. This turned out to be a red herring as the other four gamma-ray 
detections of novae all occurred in systems where the companion is a low-mass star with a puny stellar wind, yet 
the HE gamma-ray observations showed similar behaviour. The HE gamma-ray spectra are soft, peak below 1 
GeV (VHE observations have given only upper limits ]129]), compatible with either inverse Compton scattering 
or pion decay. With a luminosity « 3 x 10^® W, about 0.1% of the kinetic energy of the ejecta is channeled into 
gamma rays. Their orbital periods are not established yet ]130, 131, 132] but are likely to be of order of a fraction 
of a day, implying binary separations of a few solar radii. The envelope in these novae expands well beyond the 
binary system into the interstellar medium. However, the HE gamma-ray emission lasts « 2 weeks, leaving little 
time to sweep up material in the tenuous interstellar medium compared to supernova remnants (or symbiotics). 

Various observational features suggest the coexistence in novae of outflow components with different velocities 
]133]. “Internal” shocks between these components may be strong enough to enable particle acceleration, and dense 
enough to reprocess the X-ray thermal emission generated by shock-heating into the optical band, as in type Hn 
supernovae ]133, 134]. Higher ejecta velocities are expected from heavier white dwarfs since the critical pressure 
is reached for smaller amounts of material deposited on the surface by accretion ]125]. The companion can also 
play a role in shaping the ejecta and the formation of shocks ]135]. Perhaps, all novae produce HE gamma-rays 
and only the closest ones are detected ]136]. 

Gamma-ray observations have revealed a new facet of the nova phenomenon, one that may turn out to be 
very prohtable to study since it gives a handle on diffusive shock acceleration in an environment that differs 
from the classic laboratory of supernova remnants ]51]. Simultaneously fitting the gamma-ray emission with the 
thermal emission produced by shock-heated material, sensitive to the amount of swept-up material and to the 
shock velocity, can prove very constraining for the fraction of the particles and energy injected in the acceleration 
process [127, 134]. This hope could be dashed if the origin of the shock(s) cannot be pinpointed precisely. 
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5. What’s next ? 


The gamma-ray emission detected in binaries touches on many sub-fields of high-energy astrophysics. In terms 
of processes, it touches on the physics of magnetized, relativistic outflows in gamma-ray binaries and microquasars, 
where it informs and can be informed by our understanding of gamma-ray bursts and active galactic nuclei ; it 
touches on diffusive shock acceleration in colliding wind binaries and novae, where it is comparable to observations 
and theories of supernova remnants or cluster shocks. Gamma-ray observations of binaries may yield irreplaceable 
insights: for instance, if inverse Compton gamma rays from upscattered stellar photons do trace the free pulsar 
wind (unobservable in PWN), or when gamma rays sign non-thermal processes in CWB (if deeply-embedded radio 
synchrotron emission is hidden by free-free absorption in the stellar winds) or novae (where high ejecta densities 
may hide the thermal emission from the shock-heated material [134]). 

The variety of binaries emitting gamma rays proves that particle acceleration is much more widespread than 
would have been considered a decade ago. At present, this variety is limited to the HE gamma ray regime but there 
is no reason to expect less in the VHE regime. With a ten-fold increase in sensitivity, the Cherenkov Telescope 
Array (CTA, [137]) will enable us to see deeper below the tip of the iceberg (detecting, for instance, VHE emission 
from low-mass gamma-ray binaries, flaring microquasars, or CWB [138]). VHE information is essential to constrain 
the maximal particle energy and characterize the efficiency of the accelerator [139]. For all we know, binaries may 
even add the odd contribution to the sea of Galactic cosmic rays [140]. One may also hope that these sensitive VHE 
observations will uncover emission from unexpected sources, including binaries. Faint, hard gamma-ray sources 
could be lost in the Galactic diffuse emission in the Eermi-LAT band and pop up in a Galactic plane survey with 
CTA. For instance, HESS J0632+057 has not been detected in HE gamma rays despite being a bright VHE source 
[14]. One lesson to be learned from the present observations is that the HE and VHE regimes are rarely simple 
extrapolations of each other. 

Variability is a main feature of gamma-ray emission from binaries. All of the sources in Table 1 are variable. 
Characterizing this variability constrains the timescales of the acceleration and radiative processes that are in¬ 
volved. Significant progress can be expected with CTA [137]. Its sensitivity accross a broad band of energies, from 
a few tens of GeV to tens of TeV, will enable high quality spectral energy distributions to be reconstructed on 
much shorter timescales than possible now, allowing e.g. for exquisite phase-resolved studies of the emission from 
gamma-ray binaries. However, compared to other objects, there is no escape from the additional efforts and perse¬ 
verance required to obtain a complete observational picture, identify the different regimes of gamma-ray emission, 
accumulate high signal-to-noise spectra for each regime, establish the relationships to variability at other wave¬ 
lengths, and measure the basic binary parameters when the system is poorly known (nature of the components, 
orbital period...). 

Finally, on the longer term, the frontier of gamma-ray observations will have to move to lower energies, in 
the 1-100 MeV range (LE gamma rays) [137[. This range is undoubtedly rich in sources associated with binaries. 
Emission can be expected from the very hot plasma in accretion disk coronae, from non-thermal particle injec¬ 
tion, from radioactive decay or spallation. The current sensitivity in LE gamma rays is poor compared to the 
HE/VHE gamma-ray regimes and compared to the neighbouring hard X-ray regime, where the unprecedent sen¬ 
sitivity of NuSTAR observations is providing a wealth of discoveries and fresh insights [141]. In Table 1, LS 5039, 
LS I +61°303, Cyg X-3, Cyg X-1 are already known emitters of LE gamma rays but the quality of the spectra and 
lightcurves is far from those at other wavelengths. A ten-fold improvement in sensitivity in the LE range would 
certainly revolutionize our understanding of how particles in various objects, including binaries, cross from the 
thermal to the non-thermal pool. 
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